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We present a calculation of the process B° — ► .D ?/ within a heavy-light chi- 
ral quark model. We assume that the rf has a large gluonic component, and 
its coupling is described via the glue-glue-// effective vertex. The main con- 
tribution comes from the non-factorizable part of the effective weak Lagrangian 
at quark level. We find, within our model-dependent assumptions, a branching 
ratio B r (B° —>■ D°rj') = (1.7 - 3.3) xlO -4 , somewhat below the experimental 
upper bound 9.4 x 10~ 4 . 



I. INTRODUCTION 



After the CLEO [1] results reporting a very large branching ratio for the inclusive pro- 
duction of rj' in B — > i]'X s decays, Atwood and Soni [2] proposed an interesting explanation 
of the data. The suggested mechanism is based on the subprocess b — > sg* — > sr/'g, where 
the virtual gluon g* emerging from the standard model penguin couples to rj' via an effec- 
tive gg*r)' vertex related to the gluonic triangle anomaly. The structure of this vertex was 
re-examined in [3], where the running of the effective coupling of rj' to gluons, assumed to be 
constant in [2], is also taken into account. The possibility that the gg*rj' vertex could dan- 
gerously be affected by non-controllable non-perturbative effects was soon after discussed in 
[4]. Some further criticism can be found in [5]. 

Interestingly, the production of rj' by thermal gluon fusion at RHIC and the LHC has 
recently been discussed in [6] using again the gluon-r/ vertex idea. The authors of [7] have 
also descibed a gluon fusion process to study the inclusive B — > i]'X s and the exclusive 
B — > K^rf decays: the gluon g of the gg*rj' vertex is supposed to be emitted by the light 
quark inside the B meson, while the g* comes from the b — > s penguin. The B — > Kr)' 
decay has been further investigated in the context of perturbative QCD in [8]. To assess the 
reliability of such mechanisms, a study of the gluon content of the rj' is in order. The estimate 
found in [9] results in a gluon component as large as 26%. This seems again to favour the 
possibility that the rj' can reliably be coupled through gluons in hadronic processes. 

A different scenario proposed to explain the abundant production of rf in B — > Kr{ 
data was based on the idea of a strong intrinsic charm component in the r( [10]. Although 
appealing, the b — > ccs Cabibbo-favoured process generating an rj' via its intrinsic cc requires 
too large a charm component to explain the data. Recently the charm content of the rj' has 
been thoroughly investigated [11] and it is estimated to be too small to motivate the b — > ccs 
mechanism. 

The non-leptonic decays including an rj' in the final state, such as B° — > D°r)' and B° — > 
D*°7]', are also quite interesting on both theoretical and experimental grounds. Measuring 
the branching ratios of these processes may shed light on the nature of the rj - rj' mixing 
and of the decay constants f v , f v >. From the experimental side, at present the only known 
limits on these decays are from CLEO [12]. 



1 



The simplest possible mechanism to consider for B — > Drj' decay is naive factorization 
of quark currents as in [13]. The quark models considered there do not involve gluons. In 
contrast, in the present paper we consider the possibility that a gluon mechanism could be 
at work even in a process like B° — > D°i]', where there is no gluon arising from a penguin 
diagram (which is the case in B — > Krj'). We will assume that the main contribution to 
the decay B° — > D°r]' comes from the gluonic mechanism, as illustrated in Fig. I. The rf 
couples to the gluons as shown in Fig. 2. 

In general, non-leptonic decays are described by a quark level effective Lagrangian at 
some chosen scale. This is built up by coefficients (containing all short-distance effects above 
the chosen scale) times quark operators, typically given as the product of two currents. In 
order to get an effective Lagrangian describing the physical degrees of freedom, the mesons, 
a "bosonization procedure" is required. This typically amounts to transforming (see e.g. 
[14]) the original Lagrangian by extracting collective composite fields with boson quantum 
numbers. 

Historically, bosonizing by factorization of the currents, the so-called "vacuum saturation 
approximation" (VSA) has been a starting point for bosonizing the quark operators. The 
VSA turns out to be badly broken in fT-meson decays. However, in some classes of heavy 
meson (B) decays, factorization holds up to order l/m& [15], but not for the process we 
consider. 

A priori, it would be preferable to keep the description as model-independent as possible. 
Unfortunately, some model-dependence has to be introduced within our framework. Still the 
amplitude for B° — > D rj' presented here will be less model-dependent than the amplitude 
for the process B Krj' within the same framework. 

One should keep in mind that, when the four quark operators responsible for non-leptonic 
B decays are Fierz-transformed, they will contain products of coloured currents, such that 
a virtual gluon might be emitted from a (light) quark in the B meson due to such a current. 
Furthermore, this virtual gluon might (via the effective rj'gg* vertex) be transformed into 
an 7/ and a soft gluon. In this paper, this soft gluon is assumed to make a gluon condensate 
with another soft gluon emitted from the light quark within the D meson. 

Inclusion of such gluon condensates gave, within chiral quark models (xQM) [16,17], 



2 



significant effects for K — > 2tt amplitudes, and the chiral quark model works well for K 
decays in general [18-20]. More specific, it has been shown that the inclusion of gluon 
condensates within the xQM works well in order to understand the AI = 1/2 rule and 
the CP- violating quantity e'/e [19]. (It should also be noted that the value used for the 
gluon condensate in [19] is also in fair agreement with the result obtained within generalized 
factorization [21].) This gives a strong motivation to use models for mesons containing a 
light and a heavy quark [22-24] , and with soft gluon effects in terms of gluon condensates 
added [25]. This was already used for the process D — > K°K° in [26]. 

II. FRAMEWORK FOR THE PROCESS B° -► Lfirf 

Non-leptonic decays with AB = 1, AC = ±1 are described by an effective Lagrangian 
at quark level [27] : 

Ceff = ~^v: b v ud [CMQi + cmq 2 \ , (i) 

where <5i,2 are operators given in terms of the quark fields b, c, u, d: 

Qi = drf{\ - 75)607^(1 - 75 )m 

Q 2 = c 7 M (l- 75)^(1 - 75 K (2) 

The Wilson coefficients C\ and C 2 contain all short-distance effects down to the scale 
fi = 0(mb). As our quark model calculation of hadronic matrix elements will be based 
on the Heavy Quark Effective Theory (HQET), we should in principle also calculate the 
short-distance effects from the scale 0{rrib) down to a scale 0(A X ~ 1 GeV) and evaluate 
our matrix elements there [28] (also some new operators should then appear in the QCD 
mixing). However, in this work these complications will be neglected (more on this in 
section V). 

The operators defined in eq. (2) may be Fierz-transformed. Using the identity 

Sasti-yP = Staphs + , (3) 

the operator Q 2 takes the form : 
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FIG. 1. Bosonization of the effective Hamiltonian, T = 7 M (1 — 75) 

Q 2 = 2ct a Y(l - l 5 )udt a ^(l - l5 )b + j^Q l . (4) 

Fierzing the operator Qi gives a similar expression that describes just the process B — > Dir. 

The idea is now that the decay B° — > Z} ?/ goes in two steps. We calculate the amplitude 
for _B° — > D°g*g, where (7* is the hard virtual gluon and g is a soft gluon. Then g* — > ?/(/, 
where g' is another soft gluon. The soft gluons g and g' are then assumed to form a gluon 
condensate. For B° — >■ D°g*g we may use a generalized version [18,19,26] of the VSA, and 
obtain (see eqs. (1), (2) and (4)): 

(D°g*g\£ eff \B°) = 

^ V; b V ud {(C 1 + %) [(D°g*g\^(l- l5 )u\0)(0\d^(l- l5 )b\B°) 

+ (D°\c 1 ,(l - l5 )u\0}(g*g\dr(l - ls)b\B } 

+ 2C 2 [(gD^ct^l - l5 )u\0}(g*\dt a Y(l - l 5 )b\B°) 

+ (D^*^^ - l5 )u\0)(g\dt a r(l - l 5 )b\B )}} ■ (5) 

Within our model, the four terms in (5) give rise to four diagrams. Two of them are 
shown in Fig. 1, where we have also included the r]'g*g vertex. The other two diagrams 
are similar, except for the fact that the hard gluon is emitted from the D° meson instead 
of the B° meson. The first diagram in Fig. 1 corresponds to the factorizable part of the 
Lagrangian (proportional to C\ + C<ijN c ) and the second diagram to the non-factorizable 
part (proportional to C2). The dominant contribution comes from the non-factorizable part 
because of the large ratio of the Wilson coefficients 2C 2 /(Ci + C 2 /A^) ~ 12 at /j, — m b . It 
should also be noted that C 2 is rather stable with respect to variations in the renormalization 
scale (j,. 



v g*g 1* v l^pl^v 



FIG. 2. r\ g*g - vertex 



The form factorF / , has been examined by many groups following [2,3] (see also [29]), 
and takes the form given in Fig. 2 : 

In the Fock-Schwinger gauge one replaces e* 6 — > — ^g^-G b pi/ (0) [30], so that the r]'g*g 
vertex now takes the form: 




where q is the momentum of the virtual gluon. We will use the parametrization of Hou and 
Tseng [3] where the form factor is : 



III. ELEMENTS OF THE HEAVY LIGHT CHIRAL QUARK MODEL (HLxQM) 

Some of the matrix elements in eq. (5) are known (e.g. the ones given by the decay 
constants fp and f B ), and some are model-dependent within our approach. Our calculation 
will, to some extent, be based on what we will call a Heavy Light Chiral Quark Model 
(HL^QM) [25] (see also [26]). This is a quark loop model supposed to describe interactions 
between quarks with momenta below the chiral symmetry breaking scale A x ~ 1 GeV. The 
model is similar to the ones that have been used by several groups [22-24], but it contains 
in addition soft gluons which may form gluon condensates, similar to one version of the 
chiral quark model in the light sector [18,19,17]. The key ingredient in this kind of models 
is a phenomenological term in the Lagrangian, which couples the quark fields directly to the 
meson fields. This enables one to integrate out the quark degrees of freedom and obtain a 
Lagrangian containing only meson fields [31]. A drawback with these models is that they do 
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not incorporate any confining mechanism for quarks. However, such models, as the related 
Nambu-Jona-Lasinio models [22-24,32], describe chiral symmetry breaking reasonably well, 
and that is important for our purpose. As for the light sector [18,19,16,32], the coefficients 
of various pieces of chiral Lagrangians can be calculated. 
The Lagrangian of the HLxQM is 



£ — C-HQET + £-xQM + £-Int + ••• , 



(8) 



where 



Chqet = Q v iv ■ DQ V + 0(m Q 1 ) 



(9) 



is the Lagrangian for Heavy Quark Effective Field Theory (HQEFT). The heavy quark field 
Q v annihilates a heavy quark with velocity v and mass rriQ, and is a covariant derivative 
containing the gluon field and the photon field. 

In the pure light sector, the Lagrangian of the chiral quark model (xQM) can be written 

as : 



(10) 



where m is the constituent mass (200 - 300 MeV) of the light quarks. The field x transforms 
under the unbroken SU (3) y symmetry, so that the quark fields transforming under SU (3) l x 
SU(3) R are 



Ql = £}xl and q R = £xr , 



where £ is a 3 x 3 matrix containing the Nambu-Goldstone bosons of QCD (ir,K,r)): 



(11) 



£ = e ~ m/f , where n = —<p a (x) = 

2 v 2 



ill + ns. -4 



7T 

K- 



— 2lL _i_ _2§_ 



K+ 
K° 



(12) 



Moreover, / ~ f n — 93.2 MeV and A a , a = 1 ... 8 are the Gell-Mann matrices. The vector 
and axial vector fields are 

V M = \{id,e + £ f ^0 and A, = ~ • (13) 
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The meson-quark interaction for the heavy-light system is [22-25] : 

C-int = —GH Xf Hvf Qv + Qv H v f Xf , (14) 

where Gh is a coupling constant (in [22-24] it corresponds to a renormalization constant), 
and Hf is a field containing both the singlet (CT) and triplet (l - ) meson field: 

H} = j>(H f y>f = ((P M/ ) V - i(Pv)h 5 ) P + , (15) 

where 

P+ = (l + ffl2 (16) 

is a projection operator. The index / runs over the light quark flavours u, d, s. The field 
Ps(P M ) annihilates a heavy meson with spin-parity CT(1 _ ), with velocity v. 

The way of regularizing divergent loop integrals has to be taken as a part of the definition 
of the model. One might use some kind of cut-off regularization (sharp cut-off as in [22] 
and cut-off in the proper time formalism as in [23,24]), cutting the momenta in the loop 
between A x ~ 1 GeV and some infrared cut-off /itr- The infrared cut-off could be zero as 
in [23] or ~ m as in [22,24]. Still, dimensional regularization may be used if there are no 
mass scales bigger than A x entering the loop. Within our model we will relate all divergent 
loop integrals to some physical parameter (as for instance /„■), as was done in [18,19,17]. 
This means in particular that we will treat quadratic, linear and logarithmic divergences as 
different. The type of divergence can easily be seen by inspection. 

The calculation of soft gluon effects in terms of gluon condensates is carried out in the 
Fock-Schwinger gauge. In this gauge one can expand the gluon field as : 

A;(k) = -^G p ,(0)^\k) + ... , (17) 

where k is the gluon momentum which is put to zero after the derivation is performed. The 
fact that the soft gluon tensor is appearing explicitly in this expression makes it a powerful 
tool to derive gluon condensate contributions to many processes, as demonstrated in [30]. 
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>VG w (0)A 



kP=0 



FIG. 3. Feynman rule for a quark coupling to a soft gluon. 



Since each vertex in a Feynman diagram is followed by an integration over loop momen- 
tum, we get the Feynman rule for interactions between light quarks and soft gluons shown 
in Fig. 3. From this Feynman rule it immediately follows that the coupling of a soft gluon 
to a heavy quark is suppressed by I/tuq, since the vertex to this order is proportional to 
v fl v u G atJ-u — 0, being the heavy quark velocity. 

The diagrams contributing to the kinetic term of heavy-light mesons are shown in Fig. 4. 
By including also diagrams with the insertion of the ^(V^) fields on light quark loop lines, 
one obtains the low energy effective Lagrangian for heavy mesons coupling to the light 
mesons 7r, K, tj, . . . in the usual form [33,34] : 



'Str = "Tr 



H a (ivD-A)H a + Tr H a H b v^ ba - g A Tr H a H b ^ l5 A^ ba , (18) 



where A is the mass difference between the heavy meson(s) and the corresponding heavy 
quark. The trace is only taken over the gamma matrix indices. Note that the V M field 
and the covariant derivative (containing the photon field) can be combined into a total 
covariant derivative V^ ba = D^5 ba + iV ba - The heavy-light meson fields, P5/ and P M /, have 
been rescaled by a factor of \/Mp and \JMp*. The mass dimension of the heavy-light meson 
fields is now 3/2. Typical values of the axial coupling g A found in the literature lie between 
1/4 and 2/3. 

In order to obtain eq. (18), we have to identify the square of the coupling Gh in (14), 
multiplied by some loop integrals with the coefficients gy = 1 for the kinetic term (including 
the vector term) and g A for the axial vector term. Such loop integrals contain logarithmic 
and linear divergences plus finite terms (gluon condensate terms are finite in this case). 
Details will be given in [25]. 
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FIG. 4. Self-energy for a heavy-light meson 



The logarithmic divergent integral appearing here is the same as the one appearing in the 
pure light sector [18,19,17], and is related to /„.. The quantity 8g A = (1 — g A ) is related to 
a linear divergent integral plus finite terms (including gluon condensates). Eliminating the 
linearly divergent integral (appearing as finite in dimensional regularization) and identifying 
the logarithmic divergent integral with f n as in the pure light case, we obtain a relation for 
Gh in terms of /„., the constituent light quark mass m, 5g A , and the gluon condensate [25] : 



, 2 2m(l-\8g A ) , 96 - 3tt 

" ~~ Tt — 2 i where k = — - — 

/. 2 + ^f -2^(t G ' 2 > 768 



G u = r . — 2 , where k = — . (19) 



Thus, in the formal limit where only the logarithmic divergent integral is kept, one obtains 
the simple expression [26]: 

G H ~ ^ . (20) 

The gluon condensate and the constituent quark mass can be linked to the mass difference 
between the CT and 1~ states via the matrix element of the colour magnetic operator between 
heavy states [35] . Calculating this matrix element in our model [25] gives a gluon condensate 
range (301 MeV) 4 < (^G 2 ) < (327 MeV) 4 for a constituent light quark mass between 230 
and 270 MeV. Then we find that also numerically, the simple expression (20) is close to 
eq. (19) because the second term (~ N c m 2 ) and third term (~ (^-G 2 )) in the denominator 
of eq. (19) tend to cancel. 

When calculating the diagram for B° — > g*, the virtual gluon is not soft. (The light 
quark emerging from the weak vertex carries a momentum of order 2 GeV.) In order to 
avoid mass scales bigger than A x in the loop integral, as explained earlier, one has to split 
off this high momentum in terms of a Large Energy Effective Theory (LEET) [36-38] for 
the quark propagating between the weak vertex and the r/' vertex (indicated by two arrows 
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in Fig. 1). This is analogous to splitting off the heavy quark mass within the HQET. The 
momentum of the quark is written as p = E^n + k, where E v i is the energy of the rj' and 
n M = (1; 0,0,1) is parallel to the four momenta of the rj'. Furthermore, k is the residual 
momentum which has to be less than A x ~ 1 GeV. 

The LEET tells us that the effect of a soft gluon that couples to a fast light quark is 
suppressed by 1/E n /, since the leading order coupling is proportional to n ll n v G ailv = 0. 
Neglecting 1/mg and 1/E > corrections, soft gluons can only couple to the spectator quark 
in the B and in the D meson. 

IV. THE AMPLITUDE FOR B° -► D°r)' 

The first part of the diagram to the right in Fig. 1 describes the matrix element for the 
transition B° — > g* , where the gluon is off shell and carries a momentum q, with q 2 = rriL. 
Within our model, light vector mesons V = (p,u,(j>) might couple to light quarks like the 
field V M in eq. (10) by letting V M — > (V M — p M ) [34]. Thus, since the gluon is a vector particle, 
the transition B° — > g* is, within our model, related to the transition B° — * V = (p, u, <fi). In 
other words, writing down the loop integral for B° — > g* and for B° — > V, the integrals are 
the same, only the couplings are different. Therefore we can avoid going into details about 
the loop integral because the relevant form factor for B° — > p is known from lattice and light 
cone sum rules [39]. Using this fact will reduce the model dependence of our calculation. 
Considering a vector meson dominance (VMD) induced coupling of the p to the light quarks, 
7 M m^// p (see [40] for more details), we have in the LEET approximation: 

J"°(S° -> g* d ) = (g* d \dt a Y(l - l5)b\B°) 

^^(p\dr(l-^)b\B°), (21) 

where m p is the p mass and f p = 0.152 GeV 2 . Furthermore, a and d are colour indices from 
the quark-gluon vertex and the factor 5 ad /2 comes from the trace over the SU(3) C matrices. 
This current can now be parametrized by the form factors from the LEET result for B° — > p 
[36]: 

sad f 

J^(B° - g* d ) = —^f2E ri ,( ± [ie^v v n p el - (n ■ v)s*» + (e* ■ v)n»] . (22) 
l m p 
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A relation corresponding to (21) cannot be used for the process B — > Krj'. The other part 
of the diagram, with a light quark moving in the background field of soft gluons, contains a 
logarithmic divergence, which can be related to f n , the result of the calculation is [25,26]: 

(23) 



(gD \ct 7 „(1 - 75 )«|0) = 8m2jV G l5 



47T 

m 2 N r s 



4tt 

where v' is the velocity of the D meson. Combining eqs. (6), (22) and (23), we find for the 
non-factorizable contribution to B° — > D°rj': 



•>" S C*G H fiT D E ^X S { ±,P (n ■ v)(n ■ v>) + + . (24) 

v Qm^m^ \ n-v'J \m I N c J n 



In the calculation we have used q 2 = m 2 , in the gluon propagator. In the parenthesis 
(1 + n • v/n • v') above, the term 1 is due to the diagram to the right in Fig. 1, and the term 
n ■ v/n ■ v is due to the corresponding diagram where the r/' couples to the D meson. The 
inclusion of this diagram deserves a comment. The form factor (± for the transition B° — > p 
contains a factor \/M B [36], and it will also contain a factor l/(n • v) from the LEET loop 



integral. In this case there is also a factor \JMd coming from soft gluon emission from the 
light quark in the D meson (see eq. (23)). For the second diagram we have to interchange 
M B M Dl and v <-> v'. In total we obtain (24), where (± is the B° — > p transition form 
factor. 

The factorizable contributions are described by the diagram to the left in Fig. 1. The 
B — > Dgg* part of the diagram (connected to the 7/ vertex) is finite. The second part is 
the current (D°|c7^(l — 75)1^0), which is by definition given by the decay constant fo- The 
finite diagrams have to be calculated using the LEET, and the result of the calculation is: 



Nj " 288 mm 2 , 



] D \[M~ D n ■ v' + f B \[M~ B n ■ v 



QG 2 ). (25) 

71 



V. NUMERICAL RESULTS 



The total amplitude due to the gluon mechanism considered in this paper is: 

(DS\£e ff \B°) = -^ V ; b V ud (A F + A NF ) , (26) 
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where A F and A NF are given in eqs. (25) and (24). In the LEET limit the kinematics is 
easy: 

n-v = l, n ■ v' = M B /M D , E r( = (M§ - M 2 D )/{2M B ) . (27) 

We will use M B = 5.280 GeV, M D = 1.870 GeV, nty = 0.958 GeV , G F = 1.166 x 10~ 5 
GeV~ 2 , V* b V ud = 0.037 and f n = 93.2 MeV, which can be found in PDG [41]. The value of 
a s has to be taken at the point q 2 = m 2 ,, where a(m v 2 ) ~ 0.56 [42]. The Wilson coefficients, 
Ci and C 2 , have been calculated [27] in NLO to be Ci(m b ) = -0.184 and C 2 (m b ) = 1.078 
(where the NDR scheme and A^L. = 225 MeV have been used). 

Extrapolating C\ )2 naively down to /i — A x , we find that C 2 does not change much, 
C 2 {A X ) = 1.1695, while d is more sensitive, Ci(A x ) = -0.3695. This gives 2C 2 /(Ci + 
C 2 /N c ) ~ 98 at A x . However, this number should not be taken too seriously. Using HQEFT 
between the scales mj, and A x , there are more operators entering as in [28]. However, the 
important thing for our calculation is that the coefficient C 2 of the dominating operator Q 2 
is rather stable with respect to variations of the scale. 

We extract the numerical value for the LEET form factor (± from lattice and light-cone 
sum-rule results [39], from which we obtain the value (± ~ 0.32. 

The only model-dependent parameters entering eq. (26) are the gluon condensate and 
the constituent light quark mass m. As mentioned above, the range 230 < m < 270 MeV 
for the constituent quark mass gives the range (301 MeV) 4 < (^G 2 ) < (327 MeV) 4 for the 
gluon condensate. 

Inserting all these values in (26), and including the uncertainty in C 2 , we obtain the 
decay rate : 

M 2 — \ l - 

T(B^DS) = ^MS D \( D S\^eff\B )\ 2 = (0.7-1.4) x lO^GeV, (28) 

which gives a branching ratio 

Br(B° -> L>V) = (1-7-3.3) x 10~ 4 . (29) 

The ranges in (28) and (29) come from the variation of the constituent light quark mass 
m in the range 230 - 270 MeV, and the uncertainty in C 2 , i.e. the uncertainty within the 
model. In addition, there is the uncertainty from applying this model. 
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VI. CONCLUSIONS 



In this paper we have computed the B° — > D°r]' exclusive process coupling the rj' via 
its gluonic component and exploiting the idea of a glue-glue-?]' effective vertex, connected 
to the gluonic triangle anomaly. The calculation is performed using a heavy-light chiral 
quark model containing soft gluons which may form gluon condensates. Such a model is 
an extension of constituent quark models of the type obtained by bosonizing an underlying 
Nambu-Jona-Lasinio interaction involving heavy and light quarks. 

The results found, when compared with those in [13], where a branching ratio 
Br(B° — > D°r]') ~ 0.30 x 10~ 4 is obtained, show remarkably that the gluon mecha- 
nism could be more important than the one obtained from standard factorization of quark 
currents. It would be very interesting to check this approach directly on new data on exclu- 
sive B decays to rf mesons. At the moment only an upper bound is known for this exclusive 
branching ratio, amounting to Br(B° — > D°r]') < 9.4 x 10~ 4 [41]. 

The relevance of the gluonic cloud in the ?/, here considered from a phenomenological 
standpoint, is a crucial question within the theory. The discrepancy in mass between the r/' 
state and the octet of (n, K, rj) goldstones can directly be attributed to the gluonic content 
of r/' [43]. Moreover, the Witten-Veneziano mass formula for the rj' [44], elaborated in the 
1/N language, has very recently been derived on the lattice [45] confirming in this set up the 
profound relation between m v i and the breaking of the U(1)a symmetry due to the gluon 
anomaly. 
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